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Abstract 

First-principles calculations have been performed for the exchange interaction 
and exchange force between two magnetic Fe(lOO) thin films. It is shown that 
the magnitude of the calculated exchange force is of the order of lO^^N at 
d/a < 1 and of the order of lO^^'^N at 1 < d/a < 1.5, where d denotes the 
distance between surfaces of the two thin films and a the lattice constant of 
bulk Fe. The obtained forces are sufficiently larger than the sensitivity of 
the current atomic force microscopy of 10^^^ — 10~^^N, which suggests the 
feasibility of the exchange force microscopy. 
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I. INTRODUCTION 



Surface magnetism has attracted much attention both in fundamental and apphed re- 
searches. Experimental techniques to investigate the magnetic structure of the surface having 
been developed, mainly detect the spin polarization of electrons emitted from surface and 
the average magnetization on sample surface within a sub-micron range is usually observed. 
None of the techniques has so far probed the local magnetic structure in an atomic scale. 



provide us with a way to investigate the surface structure with atomic resolution. The STM 
probes the surface of conducting materials by measuring tunneling current between the tip 
and sample, whereas the AFM senses the atomic force between them. The AFM has a great 
advantage that it can be applied not only to conducting materials but also to insulators. 
The extension of the STM to measurements of the magnetic structure, the so-called spin- 
polarized STM (SP-STM), is currently being developed by various groups with different 
approaches. 

The magnetic force microscopy (MFM)i detects the force arising from the long-range 
interaction between magnetic dipoles of the tip and sample. The typical tip-sample separa- 
tions in the MFM are of the order of more than lOnm and the spatial resolution is of the 
order of lOnm to lOOnm. It is fairly certain that improvement of the resolution might be 
made by probing the short-range exchange force rather than the long-range magnetic dipole 



In previous papersJZfQ we pointed out the possibility of the exchange force microscopy 
(EFM) which probes the short-range exchange force between the tip and sample. Since 
microscopic understanding of the exchange interaction between the tip and sample is of 
great importance and crucial in designing the EFM, we previously studied the exchange 
interaction between ferromagnetic tip and ferromagnetic sample based on a one-dimensional 
electron gas model.&i By using a tight-binding model, Ness and Gautieri calculated the 
exchange interaction between a ferromagnetic Fe tip and magnetic Cr(OOl) and Ni(OOl) 
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surfaces. Since these calculations inevitably employ semi-phenomenological parameters, it 
is highly desirable to make quantitative evaluation of the exchange interaction and the 
exchange force for more realistic systems on the basis of first-principles electronic structure 
calculations, which is the purpose of the present paper. 

The paper is organized as follows: In Sec. II we discuss the adopted model and calculation 
method. The calculated results are presented in Sec. III. The final Sec. IV is devoted to 
conclusion and discussion. 

II. METHOD OF CALCULATION 

By modeling the tip and samples to be used in actual experiments, we adopt the two 
three- atomic-layer Fe(lOO) films which are separated by the distance ci, as shown in Fig. |l|. 
Surface atoms of the one film are assumed to be facing the hollow sites on the surface of 
the other film. The translational symmetry is preserved in the surface-normal (x) direction, 
in which a set of the two three-layer films are periodically located with the five-atomic- 
layer vacuum gap (repeated slab model). For a calculation of the electronic structure of 
the adopted films, we employ the local-spin-density approximation to the density-functional 
theory, and the linear augmented-plane-wave (LAPW) method0 with Hedin-Lundquist ex- 
change correlation.0 The potential and charge density are expanded with plane waves in 
the whole space and with spherical waves inside muffin-tin spheres, and the force is calcu- 
lated by following the prescription given by Soler and Williams.0 Because bulk Fe is in the 
ferromagnetic ground state, magnetic moments in each three-layer film is taken to be in the 
ferromagnetic alignment. As far as the relative orientation of moments in the two films are 
concerned, however, we consider the parallel (P) and anti-parallel (AP) configurations in or- 
der to calculate the exchange interaction and the exchange force between the two magnetic 
films. Our LAPW calculations are carried out by changing the film-film separation d from 
1.4Ato S.OAwith an assumption that the internal atomic coordinates of each film are rigidly 
fixed. 

III. CALCULATED RESULTS 
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A. Magnetic Moments 

We discuss first the magnetic moments in the adopted Fe thin films. Figure ^ shows the 
spin magnetic moments calculated within the muffin-tin spheres of Fe atoms in the three 
layers of the upper film, which are referred to as the xi, X2 and X3 layers (Fig. |I]). The 
film-film separation d in the abscissa is denoted relative to the lattice constant of bulk Fe (a 
= 2.83A). We should note that the separation oid/a = 0.5 corresponds just to the interlayer 
distance of bulk Fe. Due to the symmetry in our model, the magnetic moments in the x'^ 
layer of the lower film are the same as those in the x„ layer of the upper film in the P 
configuration, while in the AP configuration those in the x'^ layer has the same magnitude 
but with the opposite sign as those in the x„ layer. For both P and AP configurations, the 
magnetic moments at the outer X3 layer is enhanced to be about 2.9/1^, which is almost 
identical to the surface magnetic moment obtained by slab calculations.lli0 The magnetic 
moments at the central X2 layer is about 2.3/1^, which is close to the experimental bulk 
value and to calculated one at the central layer of the slab.0 The magnetic moments of 
the X2 and X3 layers for both P and AP configurations are almost independent of the film- 
film separation. On the contrary, magnetic moments on the Xi layer change considerably 
at (i < a, where the magnetic moments reduce significantly from the surface value to the 
bulk one as decreasing d. Nevertheless, the change in the magnitude of moments of the xi 
layer becomes insignificant for d/a > 1. We notice that the magnetic moments for the AP 
configuration are always smaller than that for the P configuration. 

B. Exchange Interaction 

The total energies, Ep and E^p, for the P and AP configurations are shown in Fig. |^ as 
a function of d/a. When the films are moved from the bulk position at d/a = 0.5, the total 
energies for both Ep and E^p increase. The exchange interaction energy, E^x, defined by 
E'ex = -E'ap — Ep is also plotted in Fig. |^. In the all region (0.5 < d/a < 1.7) investigated, we 
get the positive Ecx, which means that the P configuration is more favorable than the AP 
one. This is expected from the fact that the resulting magnetic configuration corresponds 
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to a natural stacking of ferromagnetic Fe in bulk. The exchange interaction energy has an 
RKKY-type oscillationEl with a period of about d/a = 0.7 (~ 2A). The peak values of E^x 
are 0.35eV at d/a = 0.6 and 0.19eV at d/a = 1.35. Note that the exchange interaction can 
be realized even at relatively far separation at 1.4 < d/a < 1.8 ( 4.0A< d <5.0A). 

C. Exchange Force 

Figure ^ shows the atomic forces acting on atoms at the xi, X2 and xs layers of the 
upper film. The force direction is perpendicular to the film surface due to symmetry. As in 
the case of magnetic moments (Fig. the atomic force for the Xi layer has a significant d 
dependence whereas forces for the X2 and x^ layers are almost independent of d. The total 
forces summed up these forces: E = J2nF{^n), for the P and AF configurations, are shown 
in Fig. |. The exchange force defined by = -Fap — -^p has an oscillation against d/a 
with the period of about 0.7, which is identical with the period of the exchange interaction 
(Fig. 1^). The magnitudes of the first and second peaks are 1 x 10^^ N and 2 x 10^^° N, 
respectively. In the case of d/a < 1, where the relevant d orbitals must have a large overlap, 
the exchange force may arise from the direct exchange couplings between d states of the two 
films. In the case of d/a > 1, on the other hand, the exchange coupling is expected to be 
mediated through delocalized s and p electrons. 

IV. CONCLUSION AND DISCUSSION 

In summary, we have performed the first-principles calculation of the exchange inter- 
action and the exchange force between two magnetic Fe(lOO) thin films, to show that the 
magnitude of the calculated exchange force is of the order of 10^^ at d/a < 1 while it is of 
the order of 10~^°N at 1 < d/a < 1.5. It has been shown that at > a, the perturbation 
of the approaching ferromagnetic film to magnetic moments in the other film, is negligibly 
small. The calculated magnitude of the exchange force is in fairly good agreement with our 
previous result of 10"^° ~ 10~^^iV in spite of a crudeness of the one- dimensional electron-gas 
model.ffl Our result is also comparable with the exchange force of 0.4 x lO^^iV between a Fe 
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tip and Cr(lOO) surface calculated by Ness and Gautieri with the use of the tight-binding 
model. Because the sensitivity of the current AFM is of the order of 10^^^ to our 
calculation suggests the possibility of the EFM, although it would be necessary to take into 
account an actual shape of the apex tip before we deduce definite conclusions. Ness and 
Gautier have been pointed outi that magnitudes of the exchange energy and force do not 
strongly depend on the tip morphology, although tip's shape varies the tip-sample distance 
where exchange interaction and force have the maximum values. It has been also shownl^l 
from calculations adopting sharp and blunt tips that an atomically sharp tip is not necessary 
for an actual observation with a constant-force mode, although a sharp tip is indispensable 
for a high, lateral spatial resolution. These results! justify, to some extent, our calculation 
employing a thin film as a tip; a first-principles calculation adopting more realistic tips is 
under consideration. 

In actual measurements of the surface magnetic structure by using the EFM, we would 
have the two possibilities: a constant-force mode and a constant-distance mode, just as in 
the STM.0 In the constant-force mode, we detect the variation of the tip-sample distance 
due to the force difference between the parallel and anti-parallel configurations. Such a 
experimenti has been made to detect the antiferromagnetic structure from one terrace to 
the next on Cr(lOO) surface, by using the constant-current mode in the SP-STM. On the 
contrary, in the constant-distance mode, we directly detect a change of the exchange forces 
between the parallel and anti-parallel configurations at a given tip-sample separation. If such 
a measurement is performed by changing the tip-sample distance, we would face technical 
difficulties such as the snap-in effect: when the tip-sample distance comes in the point of force 
instability where the gradient of the force exceeds the magnitude of the lever stiffness, the tip 
and surface would snap into contact.0 In order to measure the RKKY-type exchange force, 
we require an exploration of the way to avoid the snap- in effect. Recently, a continuous 
measurement of the force curve between tip and sample, as approaching towards contact 
without the snap-in effects, can be provided by means of force-controlled feedback system.lli 
This technique may possibly provide the basis of the EFM measurement. 
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FIGURES 

FIG. 1. Schematic representation of the two three-atom-layer Fe (100) films adopted in our 

calculation, surface atoms of the one film facing the hollow sites on the surface of the other film. 
Open circles represent Fe atoms and the lines denote the layer planes which are referred to as the 
Xfi and layers (n = 1 — 3) . 

FIG. 2. Magnetic moments in the muffin-tin sphere of an Fe atom in the xi (circles), X2 (tri- 
angles) and 2:3 (squares) layers as a function of the film-film separation d normalized by the lattice 
constant of bulk Fe (a = 2.83A). Solid and open marks stand for the P and the AP configurations, 
respectively. 

FIG. 3. The film-film separation dependence of total energies in the P {Ep: filled circles) and 
AF {Eaf- open circles) configurations, and of the exchange energy defined by Eg^ = Eaf — Ep 
(triangles). The reference energy for Ep and Ea.f is the total energy for the P configuration at 
d/a = 0.5. 

FIG. 4. Atomic forces, F{xn), acting on the atom of the x„ layer as a function of the film-film 
separation. Solid and open marks stand for the P and the AP configurations, respectively. 

FIG. 5. The film-film separation dependence of total forces in the P {Fp: filled circles) and AF 
(Faf: open circles) configurations, and of the exchange force defined by i^ex = Faf — Fp (triangles). 
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Fig. 1, Nakamura, PRB 
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Fig. 2, Nakamura, PRB 




Fig. 3, Nakamura, PRB 




Fig. 4, Nakamura, PRB 




Fig. 5, Nakamura, PRB 




